8213 Computer N^^rks and ISDN Systems 

.1/5, Amsterdam, NL 



• 8213 Computer Na^c 
20(1 990 )Oecember^H 




XP 000161301 
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Rate ( VBR) Sources with Application 
to Charging Principles 
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Abstract. This paper proposes a method making it possible to measure traffic from variable bit rate (VBR) sources as they occur in 
the emerging Asynchronous Transfer Mode (ATM) networks. The method is a kind of mixture of a sampling process known from 
classical communication engineering and of the Karlsson scanning method, used in the telephone networks in the Scandinavian 
countries and in Germany. The link to charging comes when we choose a charging principle where (at least part of) the charge is 
proportional to the traffic generated to the network. 
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1. Introduction 

Traffic measurement serves several purposes in a conventional telephone network: One is to make a 
general surveillance of the traffic either in order to obtain statistics for the forecasting needed for the 
future extensions of the network or in order to check the promised grade of service. Another one merely 
served the purpose of providing the user with a charge, which (s)he considers fair. (A linear function of the 
traffic generated to the network would normally be considered to have this property.) 

In the Scandinavian countries (and in Germany) the so-called Karlsson scanning principle, which was 
first described by Karlsson [1] in 1937, is used. For a telephone user the principle is the following: When a 
call is set up it will be connected to a scanning process which increments a counter with regular intervals, 
whose length only depends on the distance of the call. The scanning process used to depend on some 
hardware in the subscriber office and had nothing to do with the establishment of the individual call set-up 
processes (except that at call set-up it was ensured, that each call would be assigned to a scanning process). 
The principle is still the same though the mechanism that ensures it today is different. Using this method, 
a call whose duration is shorter than the scanning interval, will either be "lucky" in which case it 
experiences no scans at all and will be free of charge, or it will be unlucky in which case it will receive a 
scan and the call will be charged for one whole scanning interval. The number of scans is used as a 
measure for the duration of the call, and obviously the scanning method introduces an uncertainty in this 
measurement. 
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436 B. Veiro / Traffic measurement on VBR sources 

In Section 2 we shall make some general remarks on charging principles, in Sections 3 and 4 we shall 
recapture some of the main features of the Karlsson scanning method and of sampling theory respectively. 

2. General Remarks on Charging Principles 

In the "clean" version of the Karlsson scanning principle only the counts obtained during the calls are 
collected, but some network operators have, chosen a version of the method, where the counter is 
incremented as soon as the call is established in order to account for the load on the network for the 
establishment of the call. Other principles could be taken into account as well, but this paper will not enter 
the discussion of whether one or the other method is more fair to the user. Such a discussion can be found 
in [2] for instance. This paper will focus only on the part of the tariff, that has to do with the traffic 
generated to the network. 

One note should be made, however: The charging principle used in the USA for instance, where the 
subscriber gets a record of all the calls made in the period covered by the telephone bill, is impossible 
when the Karlsson method is used, because no record is kept of the fate of each individual call. The 
strength of the method is, that it is cheap to implement, but the major problem is that it relies on the fact 
that there is a relationship of confidence between the network operator and the customer. If a customer 
has any complaints it is more or less impossible both for the customer and for the administration to prove 
anything, and therefore any law-suits will be very difficult to conduct. 

3. Properties of the Karlsson Scanning Method 

The following description of the Karlsson scanning method follows [3]. Using the scanning method we 
want to measure a continuous time interval starting at / = T x and ending at t — T 2 . The number of scans 
received in the duration T=T 2 -T X is a discrete random variable Z, and we define {Z = A:} as the 
occurrence of k scans within (T u T 2 ). 

We now want to confine ourselves to constant scanning intervals, because this is the way it is actually 
implemented in practice. The fact that the method works with random scanning intervals as well at the 
cost of some accuracy, is shown in detail in [3]. If we let the scanning interval be h and the integral part of 
T x /h be \T x /h\ we obtain 

{Z = 0} = {T 2 ^ \T x /h\+h), 

{ Z = k } » { h | T x /h | +kh <T 2 ^h\ T x /h )+(/:+ \)h } (3,1) 
when k — 1, 2, 

We assume that the time interval T 2 - T x has the distribution function F(t) and that the time until the 
first scan after T x has the distribution density function v h (t). We also assume that T x is independent of 
the scanning times, which means that we can assume T x = 0 without any restrictions. Then we get the 
following probabilities: 



,{2 = 0} =/ v M (t)F(,)dt, 

, (3-2) 

p{Z = k) = f'o h (t)[F(t + kh)-F(t+{k-\)h)] dt. 
J o 

These are the general formulae applicable to any arrival process distribution and any holding time 
distribution. If we now restrict ourselves to Poisson arrivals, the time elapsed from arrival of the call until 
the next scan will be uniformly distributed in [0, h]: 

0/,(0 = 1/*, 0<(^h (3.3) 
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and we get the probabilities for general holding time distribution: 




(3.4) 



What we have done above corresponds to adding a [0, h) uniformly distributed time to the holding time 
described by F(t) and then truncating this sum to an integral time. 
By integration by parts the expectation of (3.4) becomes 



So we have the important, well-known result that the original continuous distribution and the transformed 
discrete distribution have the same mean value, which proves that the scanning principle in a sense is fair 
to the user. 

Let the time interval be of length T 2 = x + kh, k = 0, 1, 2, . . . , 0 < x < h. Then the measuring error 
becomes negative (= -x) with probability 1 - x/h and positive (= h — x) with probability x/h, and for 
the mean we have 



This is another property, that makes the Karlsson scanning fair and for the method proposed for ATM 
networks we wish to retain both the zero mean error property and (3.5). 

4. Review of Sampling Theory 

In normal textbooks in communications engineering (e.g. [7]) the proof of the sampling theorem 
depends on a formulation using the Dirac delta function (8) because it has the nice property that the 
convolution of a function /(/) with S at time t 0 gives the value /(f 0 ). If it was not for some mathematical 
problems with the interpretation of the physical implementation of sampling systems this would be very 
good because the normal idea with sampling is that we want to recover the sampled signal at the receiver. 
This is not the idea in this measurement application. First of all we only want to have a very rough picture 
of what the traffic was as long as the mean error in the measurement is kept at zero and the holding time 
distribution has the same mean value as the scanning process as was described for Karlsson scanning in 
Section 3. 

The sampling error that we shall make with the method described below will be even bigger than in 
conventional sampling systems because in normal systems we try at least to make an approximation to the 
5 function, which is as good as possible, but in ATM systems we need to take an average value over a time 
period to be specified below. If we took a true sample, we would either get the maximum bandwidth 
available for the access line or we would get zero. In principle there is nothing wrong with this, but we face 
the practical problem with the samples that we have to be able to distinguish empty time slots from 
information carrying ones and for that we need at least to have time to analyse the header, which is not 
compatible with the sampling concept. 

Therefore it would be nice if a proof of the sampling theorem that did not depend on the 8 function 
could be produced and actually such a proof has already been made in [5]. This has been translated into 
English and is provided as Appendix A as the author does not know of any sources in English providing 
such a proof. 

One more thing should be mentioned: we are not really interested in sampling in the normal sense of 
communication theory because ordinary sampling is done in order to regain the original signal (as good as 
possible with the given sampling frequency), but we are only interested in having some general statistical 
properties retained in the measurement case and only an overall picture of the signal in the charging ca.se. 




(3.5) 



>-r<-*x-a(j„ 
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5. The ATM-measurement Method Proposed 

5.L Description of Method 

The principle in the measurement method is illustrated in Fig. 1. The source starts emitting cells at time 
7^ =0. With regular time intervals the bandwidth function is "sampled" and the sampled value is digitised 
into L levels where L>2. (For L = 1 the method degenerates to ordinary Karlsson scanning.) For traffic 
measurements one would probably prefer to keep a count for each level so that after the measurement we 
are able to draw a histogram showing the bandwidth distribution. For charging purposes we only need one 
counter, which we then increment with c{k\ where c is a charging function and k is the sample number 

k e 0- 2 N )- A simple proposal for c(k) would be c(k) = k { * j + k 2 , where )G {0. 1 L } (j : = 0 

corresponds to bandwidth = 0) is a level-number and k } and k 2 are integer constants with k x = 1 and 
A- 2 = 0 as the obvious choices. Some administrations might want to have a higher price penalty for high 
bandwidth utilisation and it might be desirable to have k 2 > 0 in order to give the customer bigger 
incentive to close down the connection, if it is not actually used over a longer period of time. 

Let us now turn to the sampling method. First of all we need a measurement device that is synchronised 
to the access line in such a way that it is able to detect the individual cells and it has to know the protocol 
on the line in order to be able to detect dummy cells that do not carry any information. The measurement 
itself has to be a mean value taken over some time = T K . If we have a non-linear charging function c. we 
need to keep track of each individual source because if high bandwidth utilisation is penalised via c it has 
to be on a per source basis. But even if c(k) =j is used as a charging function as proposed above, we still 
need a counter per terminal, because the individual calls may have different destinations and will therefore 

have different scanning times. Anyway we choose the charging function c(k) = j\ where y <= (0, 1 L } 

for each sample, Ar, in the measurement, because it is the simplest one. 

The minimum time we can have for 7; is the time it takes for a cell to pass the measurement 
mechanism = 7" r . If the subscriber loop has the access rate C and the cell length is b c we get 



Taking T s = T c we would only be able to measure the rates 0 or C, but such a choice would make it 
impossible to detect the individual traffic bursts. If we make 7; too long, we will measure the mean value 
all the time, and in this case we could just have applied the ordinary Karlsson scanning as described in 
Section 3. We want 7; to be so small that it can detect a little part of a burst (whatever "little" then means 
in this connection). Burst detection (and definition) is a problem in itself and is treated in [8]. We also 
want to have a value for 7; that is independent of the policing that the sources might be subjected to 
because it is still under discussion if source policing should be performed or not with e.g. Roberts [4] 
representing the opposition and Rasmussen & Sorensen [6] the defense. A fixed time corresponding to a 
fixed number of cells for the averaging is preferable for implementation reasons. This number has to be an 
integer (to make it simple) constant times the number of levels, L, we want to distinguish. This integer 
represents (in T e units) the time over which the sample should be taken. Somewhat arbitrarily we choose 



T v = b c /C. 
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the constant to be 20 (what the real value needed should be will' have to depend on experiments with real 
sources when they will be known) 

T s = b c 2QL/C. 

We have now determined the duration of the "sample" we want to take. We now need to determine the 
time between two neighbouring samples (the h in Section 3). In [3] it is shown that the distance between 
the samples does not affect the mean value of the measurement but only the accuracy. From the sampling 
theory we know that all information outside the frequency band [~Q y Q] (where the sampling period is 
h - -n/Q) is lost, so there are limits to what we can do. The scanning interval (or sampling frequency) then 
will depend on the bandwidth function and on the application, because if we have a measurement we 
might want more detailed information than we need for charging purposes. Finally, since the charge is 
proportional to the sum of the c's, administrations will find it appropriate to let the call destination be a 
parameter that influences the sampling period. The smallest possible value for h is h - r„ because having 
a smaller h would mean that we take a new sample point before we have finished measuring the previous 
one. The final determination of the sampling frequencies needed for various destinations would probably 
be done via a field trial. But in order to maintain the idea of scanning rather than sampling we shall take 



5.2. Mathematical Model 

For the assessment of the total charge c we sample the bandwidth function B(t) in Fig. 1, where the 

kth sample has value with y6{0,l L}, where j = 0 corresponds to w < \C/L % j = 1 to' \C/L < w 

<?C/L< J=2 to \C/L^w< \C/L, ... and y = £ to w > \{2L - \)C/L. The following model is 
proposed: We regard j as a random variable whose state space is {0, 1, L }. If the values of j occur 
with the probabilities 

L 

Po^ Pi . - - • • Pl where ^ /? y =l, 

we can assume, that when we have N samples the distribution will be (L+l)-nomial ie if v n + x, 
+ • * * +x L = N, then * 0 1 

/u. - Xq1x ";„ Xl1 ps°pv - • • Pi'- 

where as usual 0! is to be interpreted as 1. This function seems to be the least restrictive we can use to 
describe the sample values of the bandwidth function. 

For this (L+ l)-nomial distribution we have the expectation vector 



From this we can calculate the expectation of the charge first in terms of charging units assuming that the 
call received N scans 

L 

M = £[0-x 0 + 1-*,+ ... +L-x L ] = NYLjPi 

y = 0 

In terms of "time-bandwidth-units" (as we have in Fig. 1) we have according to the translation rule 
given above 

The variance can also be calculated, but it does not seem so interesting in this context. 
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Fig. 2. Micro-view of bandwidth function. 



5.3. Analysis of the Model 



For the length of the call (T 2 ) we can still apply the analysis in Section 3 and hence we can conclude, 
that the number of scans N during a call, when considered as a random variable, has the correct mean 
value, viz. the same as the mean value of the distribution function describing T-,. 

^Let us now look at Fig. 1 again. If we call the bandwidth function B(t). we want to prove, that 
f 0 2 B(t)di "on the average" is the same as the value /i calculated above (in the time bandwidth unit 
version), where B(t) is to be interpreted as the outcome of an experiment so that we actually consider a 
whole family of ^-functions where our B just is a sample (in the statistical sense) from this family. 

Before we begin to assess the value of this integral, some comments on the possible members in the 
^-function family are in place. If we take a microscopic look at the ^-functions, the appearance is very 
different from the kind of picture drawn in Fig. 1, because either the source emits a cell at the moment in 
which case it is the sole user of the whole bandwidth C available, or it does not emit a cell in which case 
the bandwidth is 0 (see Fig. 2). 

This picture is disrupted by the measurement method applied so the way we regard Fig. 1 is as a picture 
taken from a long distance of a series of measurements where h = 7; has been chosen. Such a picture 
might look as in Fig. 3, and this is actually the limit to which we can go (for physical reasons) evaluating 
the integral. It is clear that the integral of such a function will exist as long as it is taken over a finite 
interval. In the expectations calculated in Section 5.2 the jc/s represent the number of occurrences of the 
event w - iC/L. The >f actors in the u- calculation represent the number of area units for each count in the 
yth column and then the (x calculated corresponds to the area under the step-function shown in Fig. 3. But 
in Fig. 3 we have the limit situation where we have done the measurement as fine as possible and with the 
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Fig. 3. Micro measurement. 
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highest possible sampling frequency and therefore for this case the jjl we got is actually the same as the 
integral (by definition) and the proof is finished. 

5.4, Example 

In the CCITT SG XVIII [9] they have recently agreed that the cell length in ATM shall be 53 bytes with 
48 bytes information and a 5-bytes header. If we take a subscriber access line with the rate C- 34 Mb/s, 
we get for the time it takes for a ceil to pass the measurement mechanism 

T. - 53- 8/34 - 10 A s = 12.5 jxs 
If we take the desired number of levels L = 8. we get for the duration of a sample 
T x = T f . • 20L = 12.5 * 20 • 8 \ls = 2.00 ms. 

Normally examples are used to verify, that the theory provided actually works when it comes to 
practical applications. As the only VBR-sources we know, with a reasonable amount of confidence are 
voice sources, these will be used. Sriram and Whitt [10] report the following values for packetized voice 
traffic: 

mean talkspurt duration: 352 ms, 

mean silence duration: 650 ms (in active mode). 

ADPCM-coding at bit rate 32 kb/s. 

The probability to hit a 32 kb/s signal on a 34 Mb/s line is only 32 - \0*/34 • 10 6 = 0 94 - 10"- but then 
the source is only active a 352/(352 + 650)th of the time, i.e. 0.35, and the total probability to hit the 
signal active is only 0.33 • 10"-\ i.e. the probability not to hit is 1 - 0.33 - 10-\ and since Sriram and Whitt 
assume geometric distribution in the active mode and exponential distribution for the pauses, we can 
assume the memoryless property, i.e. the probability not to hit the signal 160 consecutive times is 
(1 0.33*10 ) =0.95. I.e. the probability to hit at least one cell is 0.05. which does not seem 
satisfactory. Of course there is no reason to provide ATM just for voice sources, because the current circuit 
switched network does that in a cost effective manner already, so to provide a good example better data is 
needed. 



6. Conclusion 



This paper proposed a measurement method for VBR sources as they appear in ATM networks The 
method is a generalisation of the well-known Karlsson scanning principle which is used for charsin- and 
measurement purposes in the Nordic countries and Germany. A proof is given which shows ^that the 
method on the average gives the correct value. In the end an example shows that the method is not verv 
well suited to voice sources, but this is to be expected as ATM networks will not be provided just to satisfy 
the need for voice communication. Future studies of more realistic VBR sources are needed for a better 
assessment of the procedure proposed. 
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Appendix 1. Translation into English of the chapter on the sampling theorem in |5] 

It follows from the definition of an entire function, the definition of the exponential type and from 
Paley-Wieners' theorem that L 2 ([ — it, it]) is a bijection onto the set of entire functions f(z) of 
exponential type it for which 

f*~\f(*)\ 2 dx< +CO. 

(for a more detailed explanation of these concepts the reader is referred to any standard work in the theory 
of complex functions (like [Al]). The set of all these entire analytic functions is called the Paley- Wiener 
space and is denoted P, i.e. 

P= |/(z) entire | J \f(x) \ 2 dx < + co, and there is a C > 0 such that for all z e C: 

|/(r)|<Ce"--'}. 

It is obvious, that P becomes a vector space over C using pointwise addition and scalar multiplication. 
Since the restriction of /<=/> to R belongs to L 2 (R\ we have a natural inner product in P given by 

<f* 8>p-<f* g>= / + dx, /(z), g(r)e/>. (A.l) 

As the restriction of the functions P to R through the Fourier transform is a bijection onto the functions 
in L 2 ([ — it, it]) and as the Fourier transform is an isomorphism. P is a separable Hilberi space which is 
isometricly isomorphic with L 2 ([ — it, tt])_ 
If P has the representation 

/(l) -^/> )e,ad '': 

where e L-([— n, it]), it follows from Plancherel's theorem (see [Al] again) that 

II /III- / + *l/(-v)| 2 dx- / +, "|*(f)l a d'-|l*llf. (A.l) 
As the Fourier coefficients for <(»G L 2 ([ — it, tt]) with respect to the orthonormal base 

\7fT e " i '" ,,jeZ } are . /( " ) = 7ib/_ + > )e "" d '- 
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we have (where the series is convergent in the L 1 sense) 

and it follows from ParsevaFs theorem, that 
ll + llf- + £ l/(»)| 2 . 
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/I = — CO 



Hence, using (A.2) we get 

11/112- L i/(»)l 2 - 

In particular we have (/(/7)) /je ^e / 2 (Z). 

As the inverse Fourier transform of e" i/, 'J!f,_ Ml /y2TT is given by 



(A.3) 



V 7 ?* -/-I- \Z2tT* 

1 1 



r e' v 'dr 



4~ r\ iix -»dt=± 

2tt7_„ 2tt 



i(-v-/i) 



tt(a — w) 2i 



i(.v-«) W _ e - ( .v-«)Tj = sin[(x - n)tr] 



, ^ . « sin ATT 
(.v-/?)ir ~ ' it(a-- /?) 



we conclude that the orthonormal base {(l/v/2^ ) e -"" | /i e Z } by the above mentioned isomorphism is 
mapped onto the orthonormal base (sin[(r - n )<tt]/( z - n )ir | n e Z } in />. Since <KO = 
T.*J z ^ yi f(n) e ""/\/2tt by this isomorphism is mapped onto /(z), we have 

sin[(z - n )tt] 



/(0 ? ? ,_ 1)7( „ )S 



(A.4) 



where the convergence of the series is defined by the metric in P. We shall make a closer examination of 
the convergence in P. 

Let g(r) = sin -nz/<nz y zeC\(0), g(0) = 1. Then is an entire function. Due to the continuit 

there exists a C > 0 such that 

\g(s)\ <C for | a | <1 and | r| < 1. 
where z = a* + i y. Furthermore we have 
*lrl 



and 



!*(--)! < 



e' 



I.W 



TT I V 



for | a- | < 1 and j r | *s 1 



for | x | ^ 1 . 



Now let n 0 € Z be chosen so that | x - n 0 | < and let a 0 = | a - /; 0 |. If | r | ^ 1 we have the estimates 
sin[<n(z-n)\ 2 +3C 



t(;-, 7 ) 



+ 



= l l?(---»)i-= £ ig(.--,, 0 -»)i 2 

/I = — 'XL 



n = — oc 



^ ^ ,^2C 2 +-^e-E 5 ^ 

v () )' -it" „ =1 (/,- i)~ 



< + 'JZ. 



and if | v| ^ 1, we have the estimates 



sin[iT( z — n )] 
<n(z- //) 



= E ls(---/i tl -/i)| 2 <e 2 ".-' -^ + ^ L 



\ — -< /_, ' T / 
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These assessments show, that if |>-| < A:, then the sum E*.TL ^ |sin|>( z - // )]/tt(z - n ) \ 2 is uniformly 
convergent. 

According to Cau^y-Schwai^ inequality we have 



11= — OO v ' 



< E i/(«)i 2 E 



sin[*ir(r — n )] 
<n(z-n) 



< -hoc. 



hence we conclude that if (f(n)) nGX e / 2 (Z), is given, the series 



sin[ir(z — n )] 
tt( z — /I ) 



is absolutely and uniformly convergent in r, when |Im-|^A\ From this we conclude, thai the 
convergence in P means that 

is uniformly convergent in every strip of the form { z e C | |Im - 1 < A- }. Hence 

/(*)- E (-!)"/(" ) J". for.eC (A.5) 

This result shows that every function / e /» can be reconstructed from the knowledge of ihe values 
{/(/?) | /? £Z}. The series representation (A.5) is called the cardinal series for f It was introduced for the 
first time in 1915 by Whittaker [A2]. Today this is known as the sampling theorem, due to the role it plays 
in communication theory (Shannon [A3] and Shannon & Weaver [A4]). 
As a consequence of the above we have the following: 

Theorem A.l. To every complex series (a n ) n ^ y ^ l 2 (Z) there exists a uniquely determined element /e P. 
such that /(/?) = a n . This element is given by 



r( \ V° sin[ir(r - n)] V? „ sin irr 



and f(x)€: L 2 (R)nC~(R) for z = x real. 

When this is translated into the argot used in the field of information theory the theorem is formulated: 

Theorem A.2. If the spectrum <£(co) of a function of the time /(/) only contains frequencies from a finite hand 
— Q, fi], i.e. <p(w) = 0 for \ | > Q. then f(t) is completely determined by the values at the sampling 
points f(nT). n e Z, where the sampling period is T = it / <2 a /7t/ 

r/ x ^° ^vSin[fl(f-/ir)] 

/l = — -x v ' 

Apart from the terminology, the only change here is that we consider the interval Q] instead of 

[ — it, it]. It is assumed, of course, that pe L 2 ([ — Q]). 

In spite of the great usability of this result a few warnings should be given here. From a mathematical 
point of view it is absurd that a time dependent function is uniquely determined by the values at the 
sampling points. The theory only works because f{t) implicitly is assumed to be analytic, i.e. /(/) can be 
extended to an entire function. In particular /(/) has a convergent series representation 

/(')= Z c„t>\ 
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Bui then the technique becomes absurd, because there is no reason whatsoever to assume that a physical or 
technical process can be described by an analytical function. On the contrary there is every reason to 
beheve that this cannot be done by a C~-function. As long as we consider the whole real axis R the 
analytical functions, whose domain contains R, lie extremely sparse in C V (R). (Just imaoine that an 
analytical function of the kind considered is uniquely determined all over its domain, when one only knows 
the values m an arbitrary small interval ] - e, «[.) Hence there is every reason to believe that the 
assumption <p(o>) = 0 for almost all to outside {-Si, Q) is never fulfilled. 

Well, one would say, why does it work anyway? The secret is that in practice one is only able to 
measure f(nT) for a finite number of sampling points, and one puts /(/;7~) = 0 for \n\> N which 
corresponds to a consideration of the time function /(/) only in the finite time interval [-NT NT] and 
in each such compact interval Weierstrass' approximation theorem already shows, that every continuous 
function can be approximated uniformly by using analytical functions (even by polynomials') and the 
absurdity has disappeared. ■>*-■> •/• 

As it was assumed above that f(nT) = 0 for | n \ > N and implicitly that f(nT) < + oo for | n I < N 
the numbers a„ =/(«T) will define a series from L 2 {Z). Since we in any case for (c„)„ mx e L 2 (Z ) shall 
have, that c„ -* 0 for n —■ + so, and as the convergence of the series 

*~ sinffl( f -„r)] 
„ ~. " a(t-nT) 

is absolute and uniform, we only need to choose N so that £„„ > M |c„ | 2 < e\ in order to ensure that the 
error we commit by putting a„ = 0 for | „ | > N and a„ = c„ otherwise, is very small. (How this is done in 
pracuce must depend on intuition, as the c n for \n\ > /V are precisely the values that have not been 
measured.) It would be obvious to put /(/) = 0 for \i\ >N+\, so that | /(/) | < e • cl where cl is a 
constant that depends only on the series (sin[0(f - nT)\/Q{t - nT)) nsy . This is equivalent to saying that 
we only consider the values of /(/) for t e [ — NT, NT). For if /(,) . 6 for \t | > N + I, then we would 
have /(z) = 0 for all ,-sZ, since /(z) was an entire function. 

In practice one should not be confused by this. The theory of Complex Functions (and by the way also 
Functional Analysis) is a powerful tool in the solution of these problems even though nature does not 
always behave like an analytic function. 
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